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In this lecture we discuss some exact results for the low-lying spectrum of the 
Dirac operator in adjoint QCD. In particular, we find an analytical expression for 
the slope of the average spectral density. These results are obtained by means of 
a generating function which is an extcnsion of the QCD partition function with 
fermionic and bosonic ghost quarks. The low-energy limit of this generating func- 
tion is completely determined by chiral (super-)symmetries. Our results for the 
slope of the average spectral density are consistent with the results for the scalar 
susceptibility which can be obtained from the usual chiral Lagrangian. 



1 INTRODUCTION 

Both from phenomenological arguments and lattice QCD simulations we know 
that chiral symmetry in QCD is spontaneously broken by the formatiorujtf | ajjiii- 
ral condensate (this issue has been discussed in several recent reviewslil&lasla). 
However, a complete analytical understanding of the underlying mechanism of 
chiral symmetry breaking is not yet available. The situation is much better 
in Supcrsymmetric Gluodynamics. In this tl^ary-it can be shown analytically 
that the chiral condensate is non-vanishing □u'tffl. One important difference 
with QCD is that in this case the fermions are in the adjoint representation. 
In both cases the partition function can be viewed as the average of a fermion 
determinant. The chiral condensate, which is the mass derivative of the free 
energy, is thus directly related to the eigenvalues of the Dirac operator. 

Generally, the Dirac spectrum cannot be obtained analytically. However, 
because the low-energy limit of theories with Goldstone bosons and a mass gap 
is uniquely determined by the pattern of spontaneous symmetry breaking, we 
expect that we will be able to derive analytical results for the low-lying Dirac 
spectrum. This program was initiated by Leutwyler and SmilgaE3 whose work 
resulted in sum-rules for the inverse Dirac eigenvalues. A complete analytical 
understanding of the low-lying Dirac spectrum came from the realization that 
it is described by a Random Matrix Theory with th e glo bal symmetries of 
QCD, also known as chiral Random Matrix Theory EJt3. This conjecture 
has been proved analytically st arfciiig from the low-energy limit of a generating 
function for the Dirac spectrumEjO. In addition to the usual fermionic quarks, 
this QCD-like partition function contains fermionic and bosonic ghost quarks 
with a mass determined by the magnitude of the Dirac eigenvalues we are 



1 



interested in. It was understood early onEM£l that the domain of validity of 
chiral Random Matrix Theory is determined by the mass scale for which the 
kinetic term of this low-energy effective theory can be neglected, i. e. when the 
wavelength of the corresponding Goldstone modes is much larger than the size 
of the box. In this domain, the thermodynamic limit can only be taken if the 
mass is decreased such that this condition remains satisfied. However, this is 
exactly what happens for the low-lying Dirac eigenvalues which scale as the 
inverse Euchde aiL-mlii me The recent work on the generating function for the 
Dirac spectrumE^HiJO made it possible to include the effect of the kinetic term 
and to study the Dirac spectrum in the physical domain with box size much 
larger than the Compton wavelength of the Goldstone modes. This allowed 
us to extract properties of the average spectral density at a scale that remains 
fbced in the thermodjpamic limit. In particular, an analytical expression for 
its slope was found E^IIZl. 

In this lecture we are interested in the QCD Dirac spectrum for fermions 
in the adjoint representation. We will calculate the slope of the Dirac spectral 
density in two different ways. First, via the scalar susceptibility which can 
be calculated by means of the usual chiral Lagrangian, and second, via the 
valence quark mass dependence of the chiral caadensate. The first method 
was originally introduced by Smilga and Stern tB for the case of QCD with 
three or more colors and fundamental fermions. The second approach relies on 
the introduction of a generating function for the resolvent of the QCD Dirac 
operator as discussed above. We will find that both methods give the same 
results. One of the advantages of the second method, which proceeds by a 
direct calculation of the average spectral density of the Dirac operator, is its 
validity for N f = 1. 

In our convention, with an anti-Hermitian Dirac operator D, the eigenval- 
ues are given by 

D(j) k = i\ k 4>k- (1) 
The average spectral density is defined as 

KA) = (£<5(A-A fe )>, (2) 

k 

where the average is over the ensemble of spectra. There are two important ob- 
servables that are directly related to the Dirac spectrum, the chiral condensate 
and the scalar susceptibility. The chiral condensate is given by, 

E,l im li m ifff, (3) 

m-oy^oo V J i\ + m 
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and the scalar susceptibility can be written as, 



K = - lim j ( r P{X)d \ ■ (4) 
y^oo y J (i\ + m) 2 y ' 

Here, and below, the Euclidean space time volume is denoted by V. If the 
spectral density near zero can be expanded as 

p(A) = p(0) + |A|p'(0) + i/(0)A 2 ..., (5) 

the chiral condensate is given by by the Banks-Casher formula El, 

S = -^, (6) 
and the infrared singular part of the scalar susceptibility is given by 

X~^log(A/m). (7) 

The chiral condensate is obtained from the first term in (|J). This term does not 
contribute to K; the infrared singular part arises from the linear term in 
Ali higher order terms in (|J) can be neglected if the chiral limit is taken at fixed 
value of the cutoff A for the integration over A. The smallest Dirac eigenvalues 
thus provide us with important information about the vacuum properties of 
QCD. 



2 QCD WITH ADJOINT FERMIONS 

The Euclidean Dirac operator for quarks in the adjoint representation is given 
by 

D = lli {d ll 8 bc + f abc Al), (8) 

where the f a bc are the anti-symmetric structure constants of SU(N C ) and the 
7 M are the Euclidean 7-matrices. Because-t.h£|-£auge fields are real this 
Dirac operator satisfies the reality relationE9iI9'tL2l 

[iD, CK] = 0, (9) 

where C is the charge conjugation matrix (C = 7274) and K is the complex 
conjugation operator. Because 

(CKf = -1, (10) 
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all Dirac eigenvalues are doubly degenerate with eigenfujictions given by <fr and 
C K ip. The linear independence of <f> and CKtp followsE9 from properties of the 
scalar product under anti-unitary transformations, 

(CK4>, 0) = ((CK) 2 cb, CK<P)* = -(CK<P, 0), (11) 

so that (CKcj), 4>) = 0. Another consequence of (^|) is that it is always possible 
to find a basis for which thc_matrix elements of the Dirac operator are arranged 
into self-dual quaternions Ei In other words, the Dyson index of the Dirac 
operator in adjoint QCD has the vafue /3 = 4. 

Because 7274 (D + m) is anti-symmetric under transposition, the square 
root of the fermion determinant is given by its Pfafhan, and the Euclidean 
partition function for N f Majorana flavors can be written as 

Nf 

Z = J DAdetV\D + M)e~ s ™ = J DA J J[ d\ f el ^c(DS fa +M fg )X^ & s 



The mass matrix M is symmetric under transposition. In the case of one 
massless flavor this theory is Supersymmetric Gluodynamics. Its properties 
have been investigated in great detail (an excellent introductory review of this 
topic is available E3) . In particular, it has been showntl that the gluino con- 
densate is non-vanishing. Since the adjoint Dirac operator in the the field of 
an instanton has 2N C zero modes, this result cannot be understood in terms 
of explicit symmetry breaking by instantons as is the case for QCD with one 
massless fundamental flavor. A non-vanishing chiral condensate in Supersym- 
metric Gluodynamics thus suggests the existence of field configurations with 
winding number equal to 1 //V»_£lQnsiderable progress was made in this direc- 
tion in several recent articles&OSo. Some evidence ior the existence of such 
configurations has been found in lattice QCD as wellc3. 

3 DIRAC SPECTRUM 

In this section we review some properties of the QCD Dirac spectrum. Because 
of the axial Ua(1) symmetry, {D, 75} = 0, all nonzero eigenvalues occur in pairs 
± Afc with eigenfunctions given by </>& and 75 4>k ■ If 0fc ~ 75 4>k , the corresponding 
eigenvalue is necessarily zero. This happens in the field of an instanton. 

For eigenvalues much larger than Aqcd, we expect that the gauge fields 
do not significantly modify the Dirac spectrum so that its spectral density is 
given by a theory of noninteracting quarks, 

p{\) - VX 3 for A -► 00. (13) 




Figure 1: Schematic picture of the QCD Dirac spectrum. The arrow denotes the region of 
interest of the present article. 



In the case of spontaneous broken chiral symmetry, the chiral condensate is 
nonzero if the thermodynamic limit is taken before the chiral limit. This can 
be understood in terms of the existence of a tower of roughly equally spaced 
eigenvalues (indicated by the wavy curve in Fig. 1) with the smallest nonzero 
eigenvalue at about one average spacing from zero, 

' 7 d4) 



mm p(0) EV 

Such accumulation of eigenvalues near zero does not occur in the free 
theory. It is only possible if the Dirac spectrum near zero is dominated by 
the interactions of the theory. Strong interactions give rise to repulsion of the 
eigenvalues which, viewed as positions of particles, condense into a Wigner 
crystal. This phenomenqii_has been studied in great detail in the context of 
Random Matrix Theory ac3. As will be explained next, the smallest QCD 
Dirac eigeojpJues are correlated according to chiral Random Matrix Theory 
(chRMT)BE. 

-important energy scale in the Dirac spectrum is the Thouless energy 
3. This is the quark mass scale m c , for which the Compton wavelength 
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of the corresponding Goldstone boson is equal to the length of the box, i. e. 

m r T, 



F 2 



= L\ (15) 



For m <C m c the kinetic term of the Goldstone modes can be ignored. In this 
domain, all theories with the same pattern of chiral symmetry breaking and 
ss_gap are equivalent (several explicit examples have been constructed 



j). In partimlar thp Dirac spectrum below m c is given by a chiral 
Random Matrix Theory Ejo. For adjoint fermions this is a random Dirac 



operator with quaternion real matrix clcmcnts with a probability distribution 
that includes the fermion determinantJjjl .J&iSi haSj-hee-p confirmed by nu- 
merous lattice QCD simulationsBl^'y'BpBB&3M; (a complete list of 
references can be found in a recent reviewu). The basis for the predictive 
power of Random Matrix Theory is universality: the fluctuation properties of 
the eigenvalues on the scale of the average level spacing are not sensitive to a 
wide class of large modifications of the probability distribution. In chiral Ran- 
dom Matrbc Theory, this W4s_first shown for QCD with three or more colors 
and fi i n da m on t a.l fcr m i on s r^H 43 ! and only more recently for the case of adjoint 
fermions 00000. 

Should we also expect an accumulation of small Dirac eigenvalues for one 
massless Majorana flavor? Let us consider the chiral condensatc, 

W-37(?BCTsnc«»+™>)- (16 > 

For m — ► the partition function is dominated by configurations with zero 
topological charge, whereas the numerator in ( |l6|) obtains its main contribution 
from the v = 1/N C = v /N c configurations for which the sum over eigenvalues in 
( |l6| ) can be approximated by the = term. With a cancellation of the l/m 
faetor from the chiral condensate and a faetor m in the fermion determinant, 
we thus find 

<AA) = ijn^w (17) 

The average fermion determinants for v = and 9 = 1 differ by a faetor 
1/(AA)V. How can we understand this? Our explanation is that the zero 
eigenvalue repels all nonzero eigenvalues such that their average position (de- 
noted by Afc) is shifted away from zero by exactly one half average eigenvalue 
spacing A A. We thus have 
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Explicit chiral symmetry breaking by instantons also requires the accumulation 
of small nonzero eigenvalues exactly as happens m the case of spontaneous 
breaking of chiral symmetry when the chiral condensate is given by the Banks- 
Casher formula iB. These results are in agreement with the analysis based on 
finite volume partition functions for Nf — 1 for which the contributions to 
the chiral condensate from the different topological sectors can be obtained 
analyticallyE3o. 

4 CHIRAL LAGRANGIAN AND SCALAR SUSCEPTIBILITY 

The adjoint QCD partition function is invariant under 

*-(:MV0(:)- (i9 > 

where U G SU(Nf) (a Ua(1) axial symmetry is broken by the anamaly). The 
gluino condensate is a color singlet with flavor structure given byEJ 

(X f CX 9 ) = 5 f9 E, (20) 

so that the SU (N f) flavor symmetry is broken to O(Nf). The chiral La- 
grangian corresponding to this pattern of chiral symmetry breaking can be 
constructed in the Standard.. way. To lowest order in the momenta and the 
quark masses it is given byE3 

Z(M)= [ dUe~ I d 4 -x[^-Ttd ll Ud IM U- x -^'Tv(MU+M^U- 1 )} ^ 

Juesu(N f )/o(N f ) 

This partition function is invariant under the flavor transformations 

U->VUV T , M — ► V*MV\ (22) 

as required by the transformation properties of the QCD partition function. 
Such invariance arguments are very powerful and even allow us to. construct 
the low energy limit of adjoint QCD at nonzero baryon densityEH, but this 
topic will not be addressed today. 

To calculate the scalar susceptibility it is convenient to introduce scalar 
sources in the mass matrix of the QCD partition function 

M = mS f9 + s a T a , (23) 

where the T a are the symmetric generators of SU (N f). The scalar susceptibil- 
ity is then given by the second derivative of the QCD partition function with 
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respect to the scalar sources 



K ab = ld Sa d Sb \ Sa=0 ^ logZ(M), 



= --( Tr— = -S ab K. (24) 

V\ D+mD+m/ 2 V; 

Since the scalar sources and the mass matrix have the same transformation 
properties in the QCD Lagrangian, they also enter in the same way in the 
chiral Lagrangian. To lowest order in chiral perturbation theory, we expand 
U = cxp[i2Tr a T a / F] to second order in the ir ficlds, 

l(U + U- 1 ) = l--^7r k 7T l T k T l , (25) 

so that the scalar susceptibility is given by 
AY 2 1 

— — — r £Y(T a T^T^ )Tr (^pbrpmrpn ^ / 



K ab = _Tr(r a T fc r')Tr(T b r"T")<- / &x#yK k {x)n(x)* m {v)* n {v)) 
S 2 



1— loop 

- ! ( .__, j ?4 Tr(T a {T k ,T l })Tr(T b {T k ,T 1 }) log(A/m). (26) 
Carrying out the traces in flavor space we 

= -S* * 4 ( ^" 2)( ^ +4) log W A). (27) 
l28ir 2 F 4 N f 8V ' ' y ' 

The scalar susceptibility can also be calculated for QCD with fundamental 
fermions and two colors or for QCD with fundamental fcrmions and three of 
more colors with Goldstonc manifold given by SU(2N f )/Sp{2N f ) and SU(Nf), 
respectively. These three cases can be distinguished by the value of the Dyson 
index given by /3 — 4, f3 = 1 and /3 = 2, in this order. For Dyson index 
different from /3 = 4, the only change in ( plf ) is the replacement (N f ±-.4) /4 — > 
(N f + The case (3 = 2 was first analyzed by Smilga and Sternta. 

If the average spectral density has a simple linear expansion as given in (^]) , 
the slope of the average spectral density follows immediately from the result 
for the scalar susceptibility. It is given by 

(P' (0)} _^ 2 (N f -2)(N f +0) 

V ~ 16n 2 f3N f F 4 ' 1 8j 

where we have included the dependence on the Dyson index. In the next 
section we will present a derivation for the slope that does not rely on this 
assumption ([|), and moreover, is also valid for N/ = 1. 
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5 GENERATING FUNCTION FOR THE DIRAC SPECTRUM 

In order to obtain a generating function for the Dirac spectrum one has to 
extend.lhe QCD partition function with additional fermionic and bosonic ghost 
quarksLL3, 



The resolvent is then given by 



(29) 



^) = \ (E = £ O.A^Z^J). (30) 

The average spectral density follows from the discontinuity of the resolvent 
across the imaginary axis 

^l = i-(E(iA + e)-E(*A-c)) 

= J_(E(iA + e) + E(-tA + e)), (31) 

Z7T 

where the second equality follows from the relation E(z) = — E(— z). Notice 
that det -1 / 2 (D + z) cannot be written as a Gaussian integral. The minimal 
generating function thus requires the introduction of ghost determinants as 
in eq. (p9|), corresponding to one complex bosonic ghost quark and a pair of 
Majorana ghost quarks. 

In the sector of fermionic quarks, the symmetry is broken by the formation 
of a chiral condensate with flavor structure as given in (20). The pattern of 
symmetry breaking is thus given by 

SU{N f + 2) ^ 0(N f + 2). (32) 

In the sector of bosonic quarks, the quadratic form in the action can be written 
as 

i>* L \ ( m a^dA = / <f)* L \ (a^ \ ( (j> L 



2 \<72<Pr) \-ma 2 J \a 2 (/)R 

-f ^ L \ ( m(J 2\ f 4>l 

2 \024> R ) \~ma 2 J \o- 2 (l) R 
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(33) 



Here, = (l,i<Jk) and (1^ = 0^ + f a bcA a ^. The kinetic term is invariant 
under U(2) transformations. However, thc ficlds no longer occur in complex 
conjugated pairs after this transformation. But notice that the kinetic term 
is also invariant under the symmetry group U* (2, R), which does not affect 
the reality properties of the quadratic form. The chiral condensate given by 
the mass derivative of the partition function as well as the mass term in the 
partition function are only left invariant by the subgroup Sp(2). The Gold- 
stone manifold corresponding to the sector of bo sopiii quarks is thus given by 
U*(2)/Sp(2) and has only one degree of freedomEjO. 

In the sector of fermionic quarks we could extended the unitary symmetry 
to Gl(Nf + 2), but a noncompact symmetry group would lead to an effective 
partition function with an incorrect small mass expansion. The symmetry 
group of the sector of fermionic quarks should thus be U (N f + 2). In addition, 
the generating function is invariant under supersymmetry transformations mix- 
ing fermions and bosons. The full symmetry group is thus given by the graded 
Lie group Gl(Nf + 2|2) which is broken spontaneously to the ortho-symplectic 
graded Lie group Osp(2\Nf + 2). Thc Goldstone manifold is then given byEa 
the maximum Riemannian submanifold of Gl(N/ + 2\2) / O sp(2\N f + 2) with 
fermionic sector given by U(Nf+2)/0(Nf+2) and bosonic sector U*(2)/Sp(2). 

We denote this manifold by G/H. 

The low-energy effective partition function is given by 

Z ca = / ^dUe-I dixC ^, (34) 

JUeG/H 

where 

£ ff = ^Str^L^C/- 1 ) - |str(M(C/ + t/" 1 )) + mg*g + ad^ d^ , 

(35) 

and U = exp(i2$/F). The last two terms in (|35| ) represent the mass term and 
the kinetic term of the super- rf fiavor-singlet field $o = Str(<i>). This partition 
function has the same transformation properties as the generating function 
( p9| ) including the explicit breaking of an axial d?Z(l 1 1) symmetry. The mass 
matrix is given by M. = diag(m, • • ■ , m, z + J, z + J, z, z) , with N t masses equal 
to m. This partition function has both bosonic and fermionic Goldstone bosons 
with mass M vv = V2zE/F, M vs = y/ (m + z)?, /F and M ss = v / 2mS/ J F. The 
generating function (|2^) was first introduced to study the quenched approxi- 
mation in QCDE3. For symmetry class /3 = 2 a version of the effective partition 
function based on compact supergroups was first introduced by Bernard and 
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Goltermanii. It gives the correct perturbative expansion, but thc nonpertur- 
bative integrations over U are not reproduced correctly. Notice that a function 
is not determined by its asymptotic expansion. 

To lowest order in chiral perturbation theory, the resolvent is simply given 
by tadpoles coming from the differentiation with respect to the source J. Thcrc 
are three different types of contributions corresponding to three different kinds 
of mesons that can be excited by the source J. First, there are tadpoles with 
bosonic mesons of mass M vs that do not mix with the super-ry', second, there 
are tadpoles with fermionic mesons of mass M vv , and third there are tadpoles 
with bosonic mesons of mass M vv that mix with the super- 77'. We thus find 
the following result for the resolvent (pOj) 



1 



N f 



1 



(36) 



where the three terms in-between the braces correspond to the three different 
types of contributions discussed above, respectively. The trace of the propa- 
gator of the first two types of mesons is given by 



A (M 2 



-r— 



1 , r2 , M 2 

= M loe 

V ^ p 2 + M 2 16tt 2 S A 2 ' 

p 



(37) 



where A is a momentum cutoff. The propagatpx of the third type of mesons is 
more complicated but it is known analyticallyEZI. In thc limit of mo — 00, the 
trace of this propagator simplifies to 



T/ t—i 



1 p 2 + M 2 



V^[p 2 + M 2 v N f (p 2 + M 2 v ) 2 



(38) 



Using the explicit expressions for the trace of the propagators we find 
E 



E(*) = S 



1 - 



where /i = A 2 F 2 /2S. In the limit of m — > this result simplifies to 







(39) 



S(z) 



E(N f -2)(N f 



32n 2 N f F 4 



4) z 
— z log — 



(40) 
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From the discontinuity of the resolvent along the imaginary axis calculated 
using the second eguation of ([H]) we find that the spectral density in the limit 
m — > is given by 



(P(A)> _ S 



1 + 



(N f -2)(N f +m 
167T/37V/F 4 



|A| 



(41) 



V 



where the Dyson index is (3 = 4 in the case of adjoint QCD. We have also 
included the results for QCD with two colors and fundamental fcrmions (/3 = 1) 
and QCD with three or more colors and fundamental icpnions (/3 = 2). The 
latter two cases can be derived along the same lines E3tZl. This result is in 
agreement with the slope obtained from the scalar susceptibility. It shows that 
the spectral density can be expanded in powers of |A|. Finally, we wish to 
emphasize that the above derpation is also valid for N f — 1. Results from 
instanton liquid simulations are consistent with eq. (pil]). 

As an alternative to the supersymmetric method, the mass dependence 
of the resolvent can be calculated by introducing n flavors of fermionic ghost 
quarks with mass z and take the limit n — > at the end of the calculation. 
This so-called replica method was used to derive the low-energy limit pf. the 
quenched scalar susceptibility in lattice QCD with staggered fermionsc3. A 
critical comparisop_pf the supersymmetric calculation and the replica calcula- 
tion was given inEIl. However, we stress that, disregarding exceptiojaaL cases 
that the asymptotic series terminatesE3'L3, only perturbative results have 
been obtained by the replica method. 

6 CONCLUSIONS 

The QCD Dirac spectrum can be obtained from the discontinuity of the re- 
solvent of the Dirac operator. Its generating function is given by the QCD 
partition function with additional bosonic and fermionic ghost quarks. Undcr 
the assumption of maximum breaking of the axial symmetry, the low-energy 
limit of this generating function can be written down on the basis of the global 
symmetries of the theory. The leading infrared singularity of the resolvent, 
which provides us with the slope of the average spectral density, is obtained 
from a simple one-loop calculation. This result, and results for two other pat- 
terns of chiral symmetry breaking, can be summarized in a single formula that 
depends in a natural way on the Dyson index of the symmetry class. 

Our results for the spectral density are consistent with the infrared singu- 
larities of the scalar susceptibility which can be calculated by the usual chiral 
Lagrangian without relying on ghost quarks. Amazingly, the two calculations 
also agree for Nf = 1 when the scalar susceptibility cannot be calculated by 
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chiral perturbation theory. Apparently, it is possible to perform an analytical 
continuation in N f. Since, as an alternative to the supersymmetric generating 
function, the resolvent can also be calculated from an analytical calculation in 
the number of additional fermionic flavors, the agreement for Nt = 1 should 
not come as a surprise. 

Lattice QCD with two colors and staggered (fundamental) fermions is in 
the same symmetry class as QCD with adjoint fermions. Our results have 
been extended to quenched lattice QCD and an impressive agreement between 
analytical and numerical resulis. for the connected and disconnected scalar 
susceptibilities has been foundcJ. We are looking forward to a direct lattice 
calculation of the slope of the average spectral density as well. 
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